Simulation of Electron Transport through a Quantum Dot with Soft Walls 
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We numerically investigate classical and quantum transport through a soft- wall cavity with mixed 
dynamics. Remarkable differences to hard- wall quantum dots are found which are, in part, related to 
the influence of the hierarchical structure of classical phase space on features of quantum scattering 
through the device. We find narrow isolated transmission resonances which display asymmetric 
Fano line shapes. The dependence of the resonance parameters on the lead mode numbers and on 
the properties of scattering eigenstates are analyzed. Their interpretation is aided by a remarkably 
close classical-quantum correspondence. We also searched for fractal conductance fluctuations. For 
the range of wave numbers &f accessible by our simulation we can rule out their existence. 
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I. INTRODUCTION 

Theoretical investigations of ballistic transport 
through microstructures have shown that the spectral 
and transport properties of phase-coherent quantum 
systems, commonly called "billiards", depend strongly 
on the nature of the underlying classical dynamics. 1 
To date, most investigations have focused on the two 
limiting cases of systems with either purely chaotic 
or regular dynamics. However, neither of these cases 
is generic^ For the semiconductor quantum dots that 
are realized in the experiment^ a classical phase space 
structure with mixed regions of chaotic regular motion 
is expected. This is due to the fact that the boundaries 
of such devices are typically not hard walls (as in most 
theoretical investigations) but feature soft wall profiles 
for which such a "mixed" phase space is characteristic^ 
The investigation of quantum transport through soft- 
walled microstructures is the primary goal of the present 
communication. 

The mixed classical phase space results in specific trans- 
port properties. Consider, e.g., the classical escape rate 
from an open billiard. For a dot with "hard chaos" , i.e. a 
metrically transitive system, the dwell time distribution, 
or equivalently the length distribution P(l) decays 
exponentially, P(l) ~ e~ l l l (with I being the mean path 
length). In a mixed system, however, trajectories can 
be trapped in the vicinity of regular islands, leading 
to an increased length distribution which typically 
features an algebraic decay law, P(l) ~ (l/l)~ a . In those 
hard-wall billiards with shapes that allow for a mixed 
phase space, it was shown that trapped trajectories lead 
to quasi-bound states in the corresponding quantum 
transport problem and appear as isolated resonances in 
the conductancei^Si 7 . In the chaotic-to-regular crossover 
regime also so-called "Andreev-billiards" 8,9 and the 
effect of shot noise suppressioniSiiiii^ have recently 
been discussed. Furthermore, mixed classical dynamics 
was proposed as a mechanism giving rise to fractal 
conductance fluctuations (FCF)iH Several experiments 
have meanwhile been performed to test this prediction. 



First experimental data 
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appear to support 



this notion. However, in the corresponding numerical 
studies, no fractal structure in the conductance fluctu- 
ations could be found^SiLiS which, in part, has led to 
a number of theoretical works that propose alternative 
and sometimes even contradictory explanations for 
pCp i 20 J 2i 4 22 i 23 In the present paper we inquire into 
the appearance of this features for transport through 
soft-walled quantum dots. 

We calculate transport coefficients and scattering 
wave functions by solving the time-independent onc- 
particle Schrodingcr equation for a two-dimensional 
scattering device. Of particular interest is the semiclas- 
sical limit of transport, where the Fermi wavelength 
Xp = 2tt/\/2Ef (in a.u.) is much smaller than the linear 
dimension D of the quantum billiard, <C D. This is 
because the ratio (Xp/D) determines the resolution with 
which quantum mechanics can resolve the underlying 
(mixed) classical phase space. Note, however, that 
the semiclassical limit for transport in the leads of 
width d, Xp <C d, cannot be reached. The asymptotic 
incoming and outgoing scattering states thus remain 
in the quantum regime. In the limit of large Fermi 
energies Ep quantum transport simulations are quite 
demanding. In order to reach the high energy regime 
we employ the Modular Recursive Green's Function 
Method (MGRM),^^ which is a variant of the standard 
recursive Green's function approach^ 6 , suited for small 
wavelength. 

For a detailed analysis of the influence of the mixed 
classical phase space on the quantum scattering problem 
we compare the classical Poincare surface of section 
with the Husimi distribution derived from the scattering 
wavefunctions. We are thereby able to classify the 
isolated conductance resonances and find that a recently 
suggested classification^ in terms of scattering states 
corresponding to classically regular or trapped trajecto- 
ries, has to be extended to include resonances which are 
associated with unstable periodic orbits. The quantum 
counterpart of these orbits (commonly called "scars"—) 
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emerge in the wavefunction densities which we calculate 
numerically. These findings are supported by very recent 
experimental investigations of a soft-wall microwave 
billiard for which scarred wavefunction have been, 
indeed, observed^ We demonstrate that the resonances 
in conductance follow the characteristic asymmetric 
Fano lineshape^ and perform a statistical analysis of 
the distribution of resonance widths, amplitudes and 
Fano asymmetry parameters. Finally, we inquire into 
the occurrence of fractal conductance fluctuations. 
This paper is organized as follows. In Sec. [H] we 
present the scattering device investigated in this work. 
Section IIIII is dedicated to a discussion of the isolated 
conductance resonances and the corresponding wave- 
functions. The Fano profile of the resonances is analyzed 
in Sec. IIVI In Sec. [3 we discuss fractal conductance 
fluctuations and Sec. IVII finally gives a summary of the 
results. 



II. CLASSICAL DYNAMICS 

We consider in the following a stadium-shaped quan- 
tum dot with a semicircle of radius R and leads attached 
to the straight section (see Fig. QJi) . By focusing on scat- 
tering states with odd parity under reflections (y — > — y), 
the geometry can be reduced to one semicircle with at- 
tached leads with half of the original width. The resulting 
cavity boundary, including the added soft-wall potential, 
is depicted in Fig.nj),c. We inject electrons from the left 
into this system and study the transmission and reflec- 
tion probabilities classically as well as quantum mechan- 
ically. Classical simulations are performed by calculating 
many different trajectories of electrons with Fermi en- 
ergy Ep which enter the dot at x = —R. The initial 
positions across the lead width are uniformly distributed 
with an angular distribution P(6) ~ cos(0). The ballis- 
tic quantum scattering problem is solved with the Mod- 
ular Recursive Green's Function Method (MRGM)i 24 i 25 
We calculate scattering wave functions and the S matrix 
of the system at different Fermi energies Ep = kp/2 (in 
a.u.), where N — int(dkp/iT) modes are transmitting in 



the leads. The total transmission T is then given by 

jV iV N 
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FIG. 1: (a) Bunimovich stadium with leads attached at the 
straight segments, (b) Reduced scattering geometry consist- 
ing of a semi-circle of radius R — 1 which is attached to a 
rectangle of length L — 2R and width d = 0.3. The solid 
lines represent hard-wall boundary conditions and the dashed 
curves are contour lines of the soft-wall potential [for details 
see Eq. 12.31 1. (c) Illustration of the potential surface corre- 
sponding to contour plot (b). 
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(2.1) 



with t mn being the transmission amplitudes from incom- 
ing mode m to outgoing mode n. According to the Lan- 
daucr formula, the conductance is obtained as 



G = 2 x T . 



(2.2) 



Atomic units (h = |e| = m c ff = 1) will be used unless 
otherwise stated explicitly. 

If we choose a zero potential inside the structure and 
hard-wall boundary conditions, the scattering device cor- 
responds to the open Bunimovich stadium billiard which 
is prototypical for purely chaotic dynamicsi22i2i This be- 
havior changes drastically if a soft wall profile is intro- 
duced (see dashed contour lines in Fig. ^ an d Fig.^), 
given by: 



V(f) - 

G(d-y) 

+ &(y - d) 



-A, 

min[0, B(\x\ - 0.35) 2 - A], 

-A, 

B(S - 0.35) 2 - A, 



(2.3) 

\x\ < 0.35 
\x\ > 0.35 

S < 0.35 
S > 0.35 



with S = \r— do\ and do = (0, d). V(r) rises quadratically 
in the exterior region of the semicircle. The potential 
in the rectangular region below the half-circular module 
depends only on the x-coordinate. It follows the radial 
profile of V(r), however, only for potential values below 
the equipotential line of the lead such that injection and 
emission is barrier-free. Note that in contrast to the case 
where all boundaries are hard walls, the classical scatter- 
ing dynamics with an arbitrary soft-wall profile is not in- 
variant under scaling of the electron energy Ep. In order 
to approach the semiclassical limit of pathlcngth spec- 
troscopy, we choose a scaled potential with A = Aq x Ep 
and B = Bo x Ep. This results in classically scaling in- 
variant dynamics and all quantum results we obtain can 
be compared with one and the same classical phase space 
structure. 

As a first test of the effect the soft-wall potential has 
on transport, we plot in Fig. [21 the probability distri- 
bution P(l) for classical trajectories to leave the cavity 
after a length I. We find for the chaotic case with hard 
walls an exponential decay 2 ^*^ of P(l) and for the mixed 
case due to the soft walls a power-law behavior which 
is well approximated by P(l) ~ f — 1-5 over several orders 
of magnitude. The difference to the chaotic case can be 
understood as a signature of the trapping of trajectories 
in the vicinity of the hierarchical set of regular islands 
which is typical for mixed dynamics For a more de- 
tailed phase space analysis we plot the Poincare surface 
of section (PSS) of the classical dynamics in our soft- wall 
device. At each bounce of a trajectory against the hori- 
zontal lower boundary, the position x along the boundary 
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FIG. 2: (Color online) Log-log plot of the classical dwell prob- 
ability P(l) for the cavity with soft walls (red solid line) and 
hard walls (green dashed line). 



and the projection of the momentum vector in the hor- 
izontal direction p x are recorded. Fig. [3] shows the PSS 
for an ensemble of initial conditions with several islands 
of regular motion and unstable periodic orbits in an oth- 
erwise chaotic sea. Of particular interest are trajectories 
trapped in the vicinity of the islands (see e.g., Fig. |3J). 
The trapping region of phase space is shielded from the 
surrounding chaotic phase space by partial transport bar- 
riers which are formed by cantori as well as by stable and 
unstable manifolds. 4 A partial barrier can only be crossed 
through small gaps ( "turnstiles" ) where phase space vol- 
ume is exchanged;^ Trajectories that enter the cavity 
through the entrance lead first reach the chaotic part of 
phase space and may directly exit through the exit lead. 
Some trajectories do, however, cross the outermost par- 
tial barrier through a turnstile and stay trapped inside 
for a comparatively long time since the only path to the 
exit is via another or the same turnstile. Alternatively, 
these trajectories may get trapped even deeper inside the 
next layer of the hierarchical set of transport barriers. By 
contrast, the islands of regular motion which lie at the 
core of this hierarchy are invariant curves forming com- 
plete barriers and can thus not be accessed by classical 
trajectories emanating from the leads and contributing 
to transport. 

In order to determine the shape of the partial barriers and 
the location and size of the turnstiles the analytical meth- 
ods described in Ref.l34l could, in principle, be used. It is, 
however, quite demanding to construct the partial barri- 
ers explicitly, in particular in the present case of soft-wall 
cavities where analytical solutions for trajectories are not 
available. We have therefore determined the outermost 
barrier approximately by a numerical dwell length anal- 
ysis. We scan the dwell length I inside the cavity on a 
very dense array of points in the PSS. All initial condi- 
tions which correspond to a dwell length longer (shorter) 
than a typical threshold value, Iq — 15, are assumed to 
be inside (outside) the outermost partial barrier. As a 
result we obtain the approximate partial barrier line de- 
picted in the PSS of Fig. [21 as a green dashed line. The 



location of the barrier was found to be only very weakly 
dependent on the choice of ^o- 



III. ISOLATED RESONANCES 

Quantum dynamics profoundly modifies the phase flow 
in the presence of cantori in two ways: On the one hand, 
partial barriers become impenetrable when the phase 
space volume of the turnstile is smaller than that of the 
Planck cell (~ h) i.e. the size of the minimum-uncertainty 
wavepacket. Consequently phase flow is suppressed on a 
short time scale associated with classically allowed tran- 
sitions. On the other hand, quantum mechanics opens 
up the possibility of barrier penetration of both complete 
and partial barriers by tunneling. This purely quantum 
transport channel is, however, in general slow and as- 
sociated with the time scale for tunneling. Thus, the 
outermost partial barrier with turnstiles smaller than h 
divides the phase space into two distinct regions where 
chaotic and hierarchical eigenfunctions are concentrated 
on either side. 35 Since the hierarchical region couples only 
very weakly to the leads, quasi-bound states residing in 
this part of phase space give rise to sharp resonances in 
transmission^*^ Classically, an island of regular motion 
in phase space consists of a set of concentric invariant 
curves, each of them forming a complete barrier. Accord- 
ingly, long-lived quasi-bound states reside also in islands 




FIG. 3: (Color online) Poincare surface of section for the 
half-stadium with soft walls (see Fig. ^>,c). A large number 
of chaotic trajectories (black crosses), six regular closed orbits 
(b-e),(g),(h), two unstable periodic orbits (f) and one trapped 
chaotic trajectory (i) are shown. For orientation the arrows 
denote one point of the phase space representation of each 
trajectory. The green dashed curve shows an approximation 
for the outermost partial transport barrier. For details see 
text. 
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FIG. 4: (Color online) Comparison between total transmis- 
sion T through the soft-wall (left inset) and hard-wall geom- 
etry (right inset) as a function of the Fermi wavenumber &f 
(14 lead modes are open). Narrow isolated resonances on a 
smoothly varying background are clearly visible for the soft 
wall billiard (red solid line), but are absent for the case of a 
hard wall boundary (green dashed line). 



of regular motion and account for additional narrow res- 
onances in transmission. 

In Fig. 0] we compare the conductance (i.e. transmission) 
through our device as a function of the Fermi wavenum- 
ber Uf for soft-wall and hard- wall potentials (see insets). 
Note the remarkable difference between the two graphs: 
Sharp resonances are present in the case with a mixed 
phase space (see Fig. and completely absent for the 
chaotic case. The quantum counterpart to the classical 
PSS is the quantum phase space distribution. In the 
following we analyze the Husimi distribution H(x,p) at 
the resonance energy and investigate its localization^ We 
define the Husimi distribution (HD) in direct analogy to 
the classical PSS by projection of the scattering state 
onto a coherent state <^ coherent on the lower horizontal 
boundary^ The HD reads 



H(x,p) = \(d n ip(x,y)\cj)' 



coherent 



(x,P))\ 2 = 



(3.1) 



dx'd Tl ^*{x')e lkFp{x '- x) - (1/2)kF{x '- x) '' 



where d n ip(x) = n(x) • Vip(x,0) is the normal derivative 
of the scattering state on the lower boundary and n(x) 
is the vector normal to the lower cavity wall. 
In Fig. [3] we show both wavefunctions and Husimi dis- 
tributions for two resonant scattering states for which 
the quantum-classical correspondence between classical 
trajectories (a)-(b) and the scattering wavefunction (c)- 
(d) is particularly striking. Correspondingly, the HDs of 



these scattering states reside on top of the regular islands 
within which these periodic classical orbits propagate. 
This remarkable degree of quantum-classical correspon- 
dence allows a convenient characterization of resonances 
in terms of the underlying phase space structure. The 
observed resonances fall into three classes: Resonances 
that are associated with the regular (R) or with the hier- 
archical (H) regions in phase space, and finally those that 
are associated with scars (S), i.e. unstable periodic orbits 
in the chaotic sea. R and H resonances have also been 
found in hard-walled billiards with boundaries that al- 
low for a mixed phase space- while scars are well-known 
features in bound-state wavefunctions of closed metri- 
cally transitive ("hard chaos") systemsiSi We note that 
S resonances have been very recently identified in soft- 
walled microwave billiards Typical examples are shown 
in Fig. [HI Resonances (a) and (b) clearly fall into the class 
R, for which the HD resides on regular islands. State (a) 
corresponds to the center of the two prominent stable is- 
lands whereas state (b) resides near the outer border of 
the island. The corresponding resonance widths T mir- 
ror this difference: Although both being narrow and of 
the same order of magnitude, T w 10 -5 , the width de- 
creases from the state near the outer border (b) to the 
state in the center of the island (a), indicating that the 
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FIG. 5: (Color online) (a),(b) Classical periodic orbits in the 
soft-wall billiards which are closely mirrored in the density 
\ip(x, y)\ 2 of the quantum scattering wavefunctions at the res- 
onant energies: (c) wavenumber kF = 14.0592-7r/d for incom- 
ing lead mode number m = 14 and (d) Uf = 14.002l7r/d, 
with m = 12 . Figure (e) and (f ) give the corresponding phase 
space portraits: Husimi distribution (density plot, red) and 
classical PSS (islands, black). For reference we also plot the 
outermost partial barrier (green dashed line) and the most 
prominent regular island (blue solid lines). The blue boxes 
in the lower left corner of (e),(f) indicate the size of h which 
characterizes the quantum mechanical resolution of the clas- 
sical phase space. 
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FIG. 6: (Color online) Top row: Classical stable periodic orbits (a) and (b); a transiently trapped trajectory (c) and an 
unstable periodic orbit (d). Second row: Resonant wave function densities \ip(x,y)\ 2 closely mimicking these trajectories with 
wavenumbers (a) k F = 14.92 037r/d, (b) k F = 14.006l7r/d, (c) k F = 14.02 607r/d, and (d) k F = 8.1867r/d. Third row: The 
Husimi distributions (density plot, red) and the PSS of the corresponding trajectories (island and crosses, black). Bottom row: 
Transmission resonances, numerical results (green dots) and fits to the Fano resonance formula Eq. 14. H (red line). Fitted 
values for the width V and Fano asymmetry parameter q are indicated. 



probability to tunnel out of the island is higher near the 
border than from the center of the island (for quantita- 
tive details see Fig.|SJ). Note that these results are in close 
correspondence with the findings presented in Ref. l37l 
The HD of a "hierarchical" state (c) features pronounced 
intensity in the region near the cantorus corresponding 
to a classically chaotic trajectory that gets transiently 
trapped inside the partial barrier. Compared to the 
regular states the coupling of the hierarchical states to 
the leads is stronger and therefore results in a resonance 
width r w 10~ 3 which is typically two orders of magni- 
tude larger than for the regular states. State (d) corre- 
sponds to a "scarred" wavefunction whose classical ana- 
logue is an unstable periodic "bouncing-ball" orbit. Its 
resonance width V w 10 _1 is even larger than of most 
of the hierarchical states we recorded. We note, how- 
ever, that the very limited number of such scarred states 
which we could identify prevents us from drawing defi- 
nite conclusions about the generic value of their width. 
Nevertheless, it is reasonable to guess that the relation 
found for the present system Yr Th <C should 
hold in other systems as well. 

The close classical-quantum correspondence in the phase 
structure suggests that the width of the transmission res- 
onances can be estimated from the location of the HD 



relative to the cantorus. Specifically, we decompose the 
HD into one part that lies outside the cantorus occupying 
the area A\ in the PSS and the complementary area Ai 
lying inside the cantorus. 

To quantify the overlap of the HD with these areas we 
integrate the HD according to 

U = JJ dxd P H ( x 'P^ ( 3 - 2 ) 

In Fig. 0the ratio £,a 1 /(!;a 1 + £a 2 ) is plotted for each 
resonance as a function of the resonance width T for a 
large number of resonances which we analyzed. We find 
the proportionality 



il.27±0.12 



(3.3) 



Within the statistical uncertainty the deviation from lin- 
earity is most likely not significant. An approximately 
linear dependence on + S,A 2 ) could be expected 

for H resonances and for those R resonances correspond- 
ing to islands inside the cantorus. Clearly S and the 
remaining R resonances fall outside the validity of this 
estimate. 
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FIG. 7: The ratio £ai/(£ai +£,a 2 ) of the HDs integrated over 
the chaotic area A\ outside the outermost partial barrier and 
over the remaining phase space area Ai [see Eq. 13.21 1 as a 
function of the resonance width F. 



IV. FANO PROFILE 



near each resonance in the single-mode limit. This im- 
plies T offs = and A = 1 in Eq. (f^TTfl . As soon as k F 
passes the threshold k F — 2n/d for opening up additional 
transmitting modes, Fano resonances have, in general, a 
minimum different from zero, T ^ 0. 
In order to elucidate the formation of Fano resonances we 
decompose the total transmission T{k F ) in terms of its 
contributions from different modest A Fano resonance 
observed in T(k F ) appears also as a Fano resonance in 
the channel transmission probabilities T mn at an identi- 
cal position kn and width T. In general, the values for q 
and To are however different in all the to x n channels. 
Conversely, the sum of any number of Fano resonances 
with identical T, k R will again be a Fano resonance. Using 
the semiclassical connection between mode number and 
injection angle, sin# = ran / \dk F ), a naive expectation 
would be that for the billiard geometry of Fig.^ an d low 
mode numbers to, n a large fraction of transmission is 



The isolated narrow resonances in transmission T(k F ) 
have typically an asymmetric Fano line shaped il- 
lustrated in the bottom row of Fig. Fano reso- 
nances have been observed in many different fields of 
physics, including ballistic transport through quantum 
dotsi 25 i 38 i 39 ' 40 i 41 i 42 They occur when (at least) one reso- 
nant and one non-resonant pathway connecting the en- 
trance with the exit channel interfere. The specific inter- 
est in analyzing Fano profiles is driven by their high sen- 
sitivity to the details of the scattering process, in particu- 
lar the degree of coherence in transport and the presence 
of decoherent interactions with other degrees of freedom. 
In contrast to Breit-Wigner resonances, the asymmetric 
Fano resonances are not only determined by their reso- 
nance position kfj, width T and amplitude To, but also 
by the Fano asymmetry parameter q according to 



T(k t 



T 



■offs 



-To x 



(k F -k R + qT/2) 2 
{k F - k R f + (r/2) 2 



(4.1) 



where T offs is the offset value of the resonance mini- 
mum. Note that the smoothly varying background on 
top of which the resonance is situated is thus given by: 
T b g r _ T offs + Tq _ The relative amp ii tuc i e A of the res- 
onance (i.e. the difference between maximum and mini- 
mum value of the second term in Eq. H4.1|l 1 is given by 4 - 3 
A = To x (1 + q 2 ). Due to the time-reversal symme- 
try in our system (i.e. no magnetic field or decoherence 
present) the asymmetry parameter q can be treated as 
real and is a measure for the ratio between resonant and 
non-resonant transmission amplitude. 29 For \q\ — > non- 
resonant transmission dominates resulting in a symmet- 
ric dip at the resonant position. For q « 1 the peak is 
highly asymmetric. In the absence of non-resonant trans- 
mission, i.e. \q\ — > oo, the resonance shape approaches 
that of a Breit-Wigner profile. For coherent transport 
in the low-energy regime, where only one flux-carrying 
mode is open, the transmission will vary between its max- 
imum value T{k F ) = 1 (full transmission) and T(kp) = 
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FIG. 8: (Color online) (a) Wavefunction densities for selected 
incoming mode numbers m of a typical Fano resonance [kF = 
14.111l7r/d]. For increasing m the directly transmitted (non- 
resonant) part of the wavefunction, being dominant for m = 
1, becomes less pronounced, while the relative intensity of 
the n-shaped resonant part drastically increases, (b) The 
Fano parameter q shows no correlation to the geometric mean 
value of the in- and outgoing mode number, %Jm x n. (c) 
Amplitudes A mn are approximately proportional to yjm x n. 
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mediated by direct (non-resonant) transmission and only 
a small part by resonant transport which corresponds to 
transient trapping inside the structure. Low-mode num- 
bers in the leads correspond classically to small injection 
and ejection angles which, in our specific geometry, is 
equivalent to trajectories that connect the entrance and 
exit lead without exploring much of the cavity. On the 
other hand we expect resonant trapping to dominate for 
large mode numbers m,n. This feature is illustrated in 
Fig. [H^,, where we plot the density of scattering wave 
functions with the same kp but different incoming mode 
number m. The resonant part of the wave function (in 
form of the (H-shape) is strongly suppressed for the limit- 
ing case m = 1 but increasingly pronounced for growing 
rn. It is now of interest to explore the dependence of 
the two parameters, the partial amplitude A mn and the 
Fano parameter q mn on the "geometric" variable m x n. 
The lower limit y/m x n = 1 corresponds to (almost) hor- 
izontal injection and ejection angles with predominant 
non-resonant transport (see Fig. [SJi, m = 1) whereas 
increasing mean values of the mode numbers stand for 
larger angles and therefore rising dominance of the reso- 
nant pathway. While A mn is approximately proportional 
to m x n (see Fig. U) indicating that with large injec- 
tion and ejection angle the relative amplitude increases, 
we find the remarkable result that the values q mn are 
virtually uncorrelated to m x n (Fig. |Sp) , contrary to a 
naive picture. Intuitively, one would expect the asymme- 
try parameter q, being a measure for the ratio between 
resonant and non-resonant contributions to transmission, 
to systematically increase with yjm x n, but no signs of 
proportionality between y/m x n and q can be detected 
in Fig. Eb- 

The projections of the resonant states onto the lead wave- 
function are, however, not a basis invariant measure for 
transport. We therefore explore the transmission eigen- 
values, i.e. the eigenvalues of the transmission operator 
T = W with matrix elements 



M 



(n\T\m) = ^2{n\t\k){k\t^\m) 



(4.2) 



k=l 



and n,m,k < M corresponding to the number of open 
modes. The transmission eigenvalues Ai (i = 1, . . . , M) , 
which we label in ascending order (A$ < Aj+i, Vi = 
1, . . . , M — 1), provide a channel-basis invariant repre- 
sentation. Implicit in this analysis is the assumption 
that the matrix diagonalizing T is only weakly energy 
dependent across the width of the resonance. This is 
justified for narrow and non-overlapping resonances, as 
we have verified numerically in a few cases. We point 
out the similarity of this approach to the multi-channel 
quantum defect theory employed in atomic and molecu- 
lar physicsi 4 ^^ We explore now the appearance of Fano 
resonances in Ai . For this purpose we determine Ai for a 
large number of Fano resonances and extract from each 
Aj the offset T° ffs , the amplitude Ai and the eigenchannel 
Fano parameter qi. The corresponding averages over the 
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FIG. 9: (Color online) Offset of transmission eigenvalues T"* 
(green dashed line, open circles) and amplitudes Ai (red solid 
line, squares) of the Fano resonances in the eigenmode i. All 
Fano resonances in the ensemble considered were taken from 
the energy interval where 14 lead modes are open, resulting 
in imax = 14. The resonance parameters associated with the 
eigenmode i (for ordering see text) are ensemble averaged 
within one specific eigenmode i. Note that the amplitudes 
Ai feature pronounced values only for intermediate modes i 
where the resonance offset T° does not take on the "classi- 
cal" values close to or 1. 



ensemble of Fano resonances are denoted by T° Ss and Ai . 
As is well known from random matrix theory (RMT) 4 ^ 4 ^ 
the transmission eigenvalue distribution P(A) has the 
functional form P(A) = 1/(tta/A(1 — A)) , with a prepon- 
derance of eigenvalues near A = and A = 1. Specific 
features of quantum transport are engraved in the in- 
termediate values of the A-distribution. As it turns out, 
similar conclusions apply to the properties of Fano reso- 
nances. The values for T° s , plotted as a function of the 
eigenchannel number i in ascending order (Fig.[5J) directly 

mirrors the U-shape of P(A) with a clustering of T° S in 

the interval T° s < 0.1 and T° s > 0.9. Note that only a 
few eigenchannel numbers feature intermediate values of 

T i . Precisely those intermediate channels provide the 
dominant contribution to the resonance amplitudes Ai 
(Fig. EJ) . This observation suggests a simple semiclassical 
explanation: eigenvalues A close to A = are associated 
with classical channels of "pure" reflection while those 
close to A = 1 are associated with "pure" transmission. 
Classical transmission or reflection channels correspond 
to classical path bundels 5 ^ of sufficient size in phase space 
(> h, see Fig. |5J, such that a quantum wavepacket can 
be accommodated. 5152 Conversely, intermediate trans- 
mission values correspond to a highly structured area 
of phase space where the wavepacket encompasses both 
transmitting and reflecting classical paths giving rise to 
quantum indeterminism and interferences reflected in Ai . 
Turning now to , we find again that the Fano parameter 
is uncorrelated with Ai . Instead, the values for qi appear 
to be "randomly" distributed (not shown). To quantify 
their randomness, we plot the probability distribution of 
the magnitude \q\, denoted by P(|<?|) (Fig. HTJ|) . Within 
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0.5 1 1.5 2 2.5 3 3.5 

Iql 

FIG. 10: Distribution -P(|g|) of the absolute value of the q- 
parameter in Fano resonances found in the transmission eigen- 
values Xi (the same ensemble of resonances was considered as 
in Fig. |5J. Fitting the data (histograms) with a Gaussian 
distribution yields good agreement. 

the limited statistics available, \q\ appears to be a Gaus- 
sian distributed random variable. Since q is a measure 
for the ratio of the coupling to the resonant (r) scattering 
channel, relative to the non-resonant continuum (c), 



the distribution P(|g|) peaks near |(r|i|i)| « 0, corre- 
sponding to the limit of a "window" resonance. It should 
be noted that for individual resonances \q\ is not invariant 
under the transformation from the mode representation 
{Tmn) to the eigenchannel representation (Xi). This is be- 
cause |g| depends explicitly on the ratio of the amplitudes 
[Eq. H4.3JI ] and thus on the matrix elements themselves. 
Assuming that |(c|£|z)| is a smooth, weakly varying func- 
tion across the resonance, the coupling (or overlap) of 
the wavefunction of the resonance with the entrance (|i)) 
channel function can be considered to be a Gaussian ran- 
dom number. Such a hypothesis would agree with RMT 
predictions for chaotic wavefunctionsS even though there 
is no a-priori reason for the applicability of RMT to the 
present hierarchical phase space structure. 

V. CAN FRACTAL CONDUCTANCE 
FLUCTUATIONS BE OBSERVED? 

Our calculation extends to mode numbers up to M = 
40 thereby reaching a ratio of Xp/D < 0.01. Our cal- 
culation reaches therefore further into the semiclassical 
regime than previous calculations for 2D billiards. It is 
therefore tempting to probe for the occurrence of frac- 
tal conductance fluctuations (FCF) in the transmission 
probability. According to a semiclassical argument 13 
quantum dots with mixed classical dynamics can be ex- 
pected to give rise to self-similar fluctuations (over sev- 
eral orders of magnitude) to which a fractal (i.e. non- 
integer) dimension can be attributed. Previous numer- 



ical results failed to provide unambiguous evidence for 
the presence of FCF in these systems ASiLiS This is due 
to the limited values of kp that could be reached com- 
putationally. This difficulty can be circumvented by re- 
ducing the two-dimensional scattering devices to effec- 
tively one-dimensional systems, so called quantum graph 
models. Here numerical constraints are less severe thus 
allowing to explore a regime where both isolated res- 
onances and FCF simultaneously exists* Our present 
two-dimensional calculations for transport through the 
soft- wall stadium (Fig. do not show FCF, even for 
the highest kp that was accessible by our codes, i.e., 
kp 40x7r/dor Xp/D sw 0.01. We attribute the absence 
of FCF to the fact that even such kp are still too small 
to probe the turnstiles of the cantori in classical phase 
space. This statement is supported by our observation 
that the HDs enter the outermost partial barrier only at 
resonance energies, i.e. by tunneling. Transport into the 
inner regions of the hierarchy in phase space seems to 
be suppressed by the partial barrier. This exploration of 
hierarchical phase space by way of turnstiles in cantori 
is, however, crucial for the fractal fluctuations^ 
Our present negative result for soft-wall billiards even 
for moderately large kp leaves, however, the ques- 
tion open which mechanism is at work that has ap- 
parently produced signatures of FCF in several recent 
experimentsii2ii£ii§iiLi£ These experiments were in the 
regime where only a rather limited number of transmit- 
ting modes m, n — 2 — 6 is open^ Our results clearly 
show that the presence of soft walls in the experiment can 
be ruled out as the source for FCF at moderate kp. Note 
that this observation is in close correspondence to recent 
findings which suggest that "more complicated processes 
than those predicted in the semiclassical models are re- 
sponsible for the observed behavior of FCF" , 17 



VI. SUMMARY 

We have investigated the classical and quantum scat- 
tering properties for a soft-wall billiard with mixed phase 
space representing the generic device used in experimen- 
tal realizations. By analyzing the wave function proba- 
bility density and the Husimi distribution of scattering 
states we find remarkable similarities between the classi- 
cal and quantum phase space structures. This enables us 
to classify resonant scattering states associated with reg- 
ular, trapped and instable periodic classical trajectories. 
Such a mapping of resonant scattering states is mirrored 
in characteristic differences in the width of the corre- 
sponding resonances. Our investigations reveal that the 
observed resonances in all the partial transmission am- 
plitudes T mn (kp) follow the asymmetric Fano lineshape. 
The distribution of Fano asymmetry parameters q ap- 
pears to be surprisingly uncorrelated with the injection 
and ejection angles of the classical trajectories. However, 
the resonance amplitude is approximately proportional to 
the geometric mean of the lead mode numbers y/m x n. 
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Studying the transmission eigenvalues A^ of f — ttf , we 
find that Fano resonances in A^ feature g-parameters fol- 
lowing a Gaussian distribution and amplitudes Ai that 
have substantial contributions only in the non-classical 
transmission eigenchannels^i For numerically accessible 
wavenumbers with Xf/D « 0.01 fractal conduction fluc- 
tuations (FCF) could not be detected. 
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